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Abstract. - We analyze the effect of frequency on the width of a single finger displacing a 
viscoelastic fluid. We derive a generalized Darcy's law in the frequency domain for a linear 
viscoelastic fluid flowing in a Hele Shaw cell. This leads to an analytic expression for the 
dynamic permeability that has maxima which are several orders of magnitude larger than the 
static permeability. We then follow an argument of de Gennes to obtain the smallest 
possible finger width when viscoelasticity is important. Using this, and a conservation law, we 
obtain a lowest bound for the width of a single finger displacing a viscoelastic fluid. Our results 
indicate that when a small amplitude signal of the frequency that maximizes the permeability is 
overimposed to a constant pressure drop, gigantic variations are obtained for the flnger width. 



Many problems associated with oil, plastic, and chemical industries involve the response of 
a fluid in a porous medium to a frequency dependent pressure drop. Laboratory experiments 
have shown for example that ultrasonic vibrations can considerably increase the rate of flow 
of a liquid through a porous medium . A useful way to describe such frequency-dependent 
processes is through a dynamic permeability function Expressions for the dynamic per- 
meability of Newtonian fluids were obtained long ago for deformable and isotropic porous 
media Q . More recently, an enhancement of several orders of magnitude in the dynamic per- 
meability was found for the case of Maxwellian fluids . Many other problems in the same 
industries, involve the displacement of a viscoelastic fluid by a second fluid of lower viscosity. 
The dramatic enhancement found for the dynamic permeability of a Maxwell fluid Q flowing 
in a tube, immediately suggests that similar dramatic effects might happen in problems of 
flow involving more than one fluid. 

The viscous flngering problem has been historically very important in the area of morphol- 
ogy of interfaces out of equilibrium |^ and has been a model system to describe displacement 
of viscous fluids in porous media. The viscous flngering problem is the prediction of the shape 
of the fluid interface in a two-phase flow conflned in a Hele-Shaw cell which consists of a 
pair of glass plates parallel to each other separated by a small gap. A viscous fluid occupies the 
space between the plates and it is pushed by a second fluid whose viscosity is relatively low. 
When the fluid is pushed laterally, the experiment is said to take place in a linear cell. When 
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the fluid is injected through a hole made in one of the plates forming the cell, the experiment 
is said to take place in a radial cell. Both the linear and the radial geometry have been studied 
extensively. The interface between the fluids is unstable and the structures that are formed 
when the interface destabilizes are called fingers. This is the so called Saffman- Taylor insta- 
bihty 0. Recently, viscous fingering experiments have been made with non-Newtonian fluids. 
Some of these experiments have used clays polymer solutions and lyotropic lamellar 

phases pT[| . Non-Newtonian fluids differ widely in their physical properties, with different flu- 
ids exhibiting a range of different properties, from plasticity and elasticity to shear thickening 
and shear thinning. Several studies have been made, both theoretically and experimentally in 
order to differentiate between the effects caused by different properties ]To|,[T^-p^. 

In this letter we analyze the effect of frequency on the width of a single finger displacing 
a viscoelastic fluid. To start, we derive a generalized Darcy's law in the frequency domain for 
a linear viscoelastic fluid flowing in a Hele Shaw cell. This leads to an analytic expression for 
the dynamic permeability that in agreement with results in other geometries Q has maxima 
which are several orders of magnitude larger than the static permeability. We then follow an 
argument of de Gennes to obtain the smallest possible finger width when viscoelasticity 
is important. Using this, and a conservation law, we obtain a lowest bound for the width 
of a single finger displacing a viscoelastic fiuid. What we find is that, when there is a small 
amplitude signal of the frequency that maximizes the permeability, overimposed to a constant 
pressure drop, gigantic variations are obtained for the finger width, that is, tiny perturbations 
in the pressure drop produce dramatic effects on the amplitude of the finger width. 

We start our study taking a Maxwell fiuid, which is the simplest model of a linear vis- 
coelastic fluid, and linearize the equation governing the flow. In the frequency domain this 
equation is 

— p (trU)'^ + iw) V — ri\/^v = — (1 — iujtr) \7p . (1) 

Here both the velocity v and the pressure p are in the frequency domain, that is, they are 
functions of space and frequency, tr-, rj and p are respectively the relaxation time of the 
Maxwell model and the viscosity and the density of the fluid. We solve (|^) for a homogeneous 
fluid flow in the x direction, confined between parallel plates aX z — ±1 subject to the boundary 
conditions Wa;(±Z) = 0. We obtain the velocity profile between the plates. In order to obtain 
a generalized Darcy's law, we average over the z— direction and obtain for the average flow 

/ tmi^l\ {l~ iujtr) dp 

^''^ = V-^^)^^Tx ■ 

Here = ^ jj^ ^j^g Yvmit w — > we recover the steady state Darcy's law and for 

fy, ^ we obtain the newtonian fluid case. The dynamic permeability of a single Hele-Shaw 
cell is then given by 

Figure |^ shows the real part of the normalized permeability K{u})/K{{)) versus the fre- 
quency for both a Maxwell fluid and a Newtonian fluid. The figure shows that for a Newtonian 
fluid the permeability is a monotonically decreasing function of frequency. For a Maxwell fluid 
on the other hand, there arc frequencies in which there are resonances which give peaks for 
the permeability. For typical values of density, viscosity, plate spacing and relaxation time in 
the common literature of Hele-Shaw problems, we find that this permeability can be two or 
three orders of magnitude larger than the Newtonian permeability, this is because the behav- 
ior of the flow at such frequencies is dominated by the elastic properties of the fluid. What 
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Fig. 1 - Dimensionless Dynamic Permeability vs Frequency. The continuous line is for a Viscoelastic 
Fluid with a relaxation time tr — 6s, the dotted line for a Newtonian fluid. For both lines the viscosity 
•q = 0.7p, the density p — Ig/cm? and the spacing between the plates L — 2mm 



this result means is that if we displace the viscous fluid at the frequency that maximizes the 
permeability, the fluid will flow with the least possible resistance. We find that the frequency 
which gives the largest permeability is of the order of ten Hertz. 

Figure ^ shows how the maximum of the real part of the permeability shifts towards 
higher frequencies when the viscosity is increased. So, the more viscous the fluid, the higher 
the frequency of the necessary pumping to minimize the resistance to flow. 

We now analyze the problem of a low viscosity fluid displacing a high viscosity viscoelastic 




Fig. 2 - Effect of viscosity on the dimensionless dynamic permeability for a viscoelastic fluid. For the 
dotted line the viscosity rj — 0.6p, for the continuous line the viscosity i] — Ip. For both lines the 
density p = Ig/cm^ , the spacing between the plates L = 1mm and the relaxation time tr — 0.6a 
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fluid. In particular, we analyze the case of a single finger with a time dependent velocity U(t) 
propagating into the viscous fluid. We call X{t) the finger width. It is worth emphasizing that 
we are not considering a steady state. Both U and A depend on time 



The quantity U/X gives a characteristic frequency. The viscoelastic fiuid has a character- 
istic time tr which in the case of a Maxwell fiuid is given by tr = rj/G, G being the rigidity 
modulus of the viscous fiuid. When U/X > 1/tr the viscoelastic fluid behaves like a solid 
and there is no Saffman- Taylor instability. The allowed wavelengths should all correspond to 
U/X < 1/tr [@- Therefore, the smallest possible finger width should be such that 

M = l . (4) 

Conservation of matter implies that 

U{t)Xit) = Voo{t) . (5) 

From equations (|^) and ^ we can relate the finger width X{t) and the tip velocity U(t) with 
the velocity at the extreme of the cell Voo (t) as 

XHt) = ^V^it) . (6) 

and 

U'{t)^^Voo{t) . (7) 

Experimentally, the parameter that can be controlled is the pressure difference at the 
extremes of the cell. So, given Ap(t), we can make the following calculations: 

Ap{t) ^ Ap(a;), 

Fourier Transform 



Darcy's law in frequency domain 
Inverse F.T. 



Voo{t) ^ X{t),U{t) 

In order for fingers to exist, the pressure difference should be at any moment positive. 
Simple oscillatory signals are not possible since the instability will exist for only half of the 
period. So we are thinking, for example, in signals that are overimposed to pressure drops 
that are large enough to destabilize the interface. 

We consider the simple case of an oscillatory finger. Suppose we impose a pressure differ- 
ence of the form 

Apit)^Po+Pae-'^"\ (8) 



(^)We have recently shown that the single finger solution exists for this case. Whether or not the solution is 
stable has not been shown. 
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that is, a constant pressure drop plus an oscillatory signal of frequency wq. We obtain the 
following expressions for the different steps presented above: 

Ap{lo) — V2Trpo5{uj) + V2Trpa6{u! — ujo) (9) 



K{uj) [po5{uj) + Pa5{uJ - Wo)] 



Voo{t) - ^-K{Q) + ^if(c^o)e-*-°* 



Now, for the steady state we know that 



and then we can express all of our results as 



Pa KjuJo) 

Po K{0) ' 



-iujQt 



(10) 
(11) 

(12) 
(13) 



u\t) = ul 



Pa Kjup) 

Po K{0) ' 

Pa K{uJo) 



-iujQt 



-iujQt 



(14) 



(15) 



Po K{0) ' 

Here A^s and Vss are defined in terms of trough (^ and (|^). What it is interesting to 
note, is that the behavior is totally different if the displaced viscous fluid is Newtonian or 
viscoelastic. For Newtonian fluids, the ratio is always small or equal to one and the 

variations in the width of the finger are small. In the other hand, for viscoelastic fluids, the 
can be several orders of magnitude larger than in the Newtonian case. Hence, if 



K{0) 



ratio 

the imposed signal has a frequency uq such that it maximizes the dynamic permeability, and 
therefore the ratio ^'(^^ , a small amplitude in the oscillatory part of the pressure signal, will 
result in gigantic variations in the amplitude of both the finger velocity and the finger width. 

A word of caution is needed, since when one does not consider surface tension the finger 
width and the velocity at the finger tip are not independent. It is worth noticing that since 
the surface tension has a stabilizing effect for small wavelengths any consideration of surface 
tension would give fingers wider or equal to the lowest bound that we are dealing with in the 
present work. 

We are not solving the selection problem, in order to do it and to compute finger shapes, 
the problem should be solved as in references |Q, ^ . Nevertheless the main result of the 
present paper should remain. That is, small variations in the pressure drop at a frequency 
that maximizes the permeability, will result in gigantic variations of the finger width. 
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